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ABSTRACT 

The mysteriously low CMB power on the largest scales might point to a Universe which 
consists of a multi-connected space. In addition to a suppression of large-scale power, 
a multi-connected space can be revealed by its circles-in-the-sky signature. In this 
paper, a detailed search for this signature is carried out for those three homogeneous 
multi-connected spherical space forms that lead to the smallest large-scale power. A 
simultaneous search for all occurring paired circles is made using filtered CMB sky 
maps which enhance the ordinary Sachs- Wolfe contribution. A marginal hint is found 
for the right-handed Poincare dodecahedron at f2 to t — 1-015 and for the right-handed 
binary tetrahedral space at O tot ~ 1.068. However, due to the complicated noise 
and foreground structure of the available microwave sky maps, we cannot draw hrm 
conclusions from our findings. 
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1 INTRODUCTION 

Current observational evidence strongly corroborates the 
standard big-bang model based on a Friedmann-Lemartre 
universe possessing a space-time structure R+ x M. Here 
M+ describes cosmic time, and M the three-dimensional co- 
moving space section of constant curvature K = +1,0, —1. 
It is well-known that the Einstein field equations describe 
the local space-time geometry of the universe, once the spa- 
tial curvature K and the various mass/energy densities are 
given, but do not specify its global spatial geometry, i. e. its 
topology. Thus, in cosmology we are faced with two very 
fundamental questions: 

i) What is the value of K, i. e. the curvature of the uni- 
verse? 

ii) What is the topology of the three-space Ml 

While the curvature does constrain the possible topologies, 
it does not fix it. Only if it is assumed that the space of 
the universe is simply-connected, the possible homogeneous 
three-spaces A4 of constant curvature are given by the three- 
sphere S 3 {K — +1), Euclidean three-space £ 3 (K = 0), or 
hyperbolic three-space Ti. 3 (K — —1). At present, there is, 
however, no compelling theoretical argument why the spatial 
geometry of the universe should be simply-connected, and 
thus we are confronted with the problem of "cosmic to pol- 
ogy" fsee lLachieze-Rev fc Lumined (119951) : iLevinl (I2002T) for 
a review) . 

The discovery of the temperature fluctuations ST of the 
cosmic microwave background (CMB) radiation by COBE 



in 1992 llSmoot et al.lll99S ) and the detailed measurements 
by WMAP jBennett et al1l2003fi and by other groups offer 
an ideal testing ground for a study of cosmic topology. The 
first year WMAP data show a slight tendency to a positive 
curvature of the un iverse, at least within the la-band. In 
lAurich et alJ ll2005al ) we have therefore carried out a system- 
atic comparison of the WMAP data with the predictions ob- 
tained from a representative sample of homog eneous spher- 
ical sp ace forms. As a result, we have found iAurich et alJ 
2005a) that only three multi-connected spherical spaces are 
compatible with the observed CMB anisotropy. The most re- 
markable signature of these spherical universes is a suppres- 
sion of the CMB power spectrum at large scales, in particu- 
lar a large suppression of the CMB quadrupole and octopole, 
and of the t emperature two-point correlation functio n on 
large angles llLuminet et al J 120031 : lAurich et aljl2005bl lah in 
accordance with the COBE and WMAP measurements. It 
is worthwhile to emphasize that the suppression on large 
scales is a prediction of these universes and not put in by 
hand. 



Multi-connected spaces can re veal themselves by th e so- 
called circles-in-the-sky signature llCornish et aljll998l) . In 
this paper, we carry out a thorough search of these paired 
circles for the above mentioned three spherical spaces. 
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2 SPHERICAL FRIEDMANN-LEMAITRE 
UNIVERSES 

In our search for the circles-in-the-sky signature in CMB 
sky maps, we are only interested in large-scale fluctuations 
of the CMB above horizon at recombination since these 
fluctuations should contain the fingerprint of a non-trivial 
topology. The observed relative temperature fluctuations 
5T/T = O(10 -5 ) are very small, and thus it is sufficient to 
use linear perturbation theory. The four-dimensional space- 
time metric with scalar perturbations takes then the simple 
form (in conformal Newtonian gauge) 



ds 



a{r)) [ (I + 2$)d»7 2 - (f - 2§)\dx\ 2 ] 



(1) 



where r\ G R+ denotes the conformal time and a(rf) the 
cosmic scale factor. \dx\ is the spatial distance on a given 
spherical space form A4. For vanishing <£>, $ = 0, the metric 
Q is the well-known Friedmann-Robertson-Walker metric. 
In writing the metric JQ with a non- vanishing scalar pertur- 
bation $ 7^ 0, we have assumed that the energy-momentum 
tensor T M „ satisfies Ty = for i ^ j (i,j = I, 2, 3) and that 
vector and tensor modes can be neglected in our study of 
the CMB. 

In the following, we shall assume that the total energy 
density of the universe is given by a sum of 4 contributions: 
radiation, baryonic matter (bar), cold dark matter (cdm), 
and dark energy, where the latter is identified with a posi- 
tive cosmological constant A. Then the scale factor a(rf) is 
uniquely determined by the corresponding Friedmann equa- 
tion. 

Let us now discuss the geometry and topology of the 
spherical space forms. There are infinitely many spherical 
spaces which were cl assified bv iThrelfall fc Seifertl Jl93l 
1 19321) (see also IWolj dl974l) and iGausmann et alJ (1200 lTl ) 
and are given by the quotient M — S 3 /F of the three- 
sphere S 3 under the action of a discrete fixed-point free 
subgroup r C SO(4) of the isometries of S 3 . To define 
the subgroup F, one makes use of the fact that the unit 
3-sphere <S 3 can be identified with the multiplicative, but 
non-commutative group of unit quaternions {q}. The lat- 
ter are defined by q :— w + xi + j/j + zij, (w, x, y, z) G M 4 , 
having unit norm, \q\ 2 = w 2 + x 2 + y 2 + z 2 = 1. Here, 
the 4 basic quaternions {l,i, j,ij} satisfy the multiplication 
rules i 2 = j 2 = —1 and ij = — ji plus the property that 
i and j commute with every real number. The distance 
d(qi,q?) between two points gi and g2 on S 3 is given by 
cosd(gi, q?) = wiW2 + xix 2 + yi J/2 + ziz 2 . 

The group SO(4) is isomorphic to S 3 x 5 3 /{±(l, 1)}, 
and thus the subgroup F can act on the first S 3 factor, on 
the second S 3 factor or on both factors. Since we are only 
interested in orientable homogeneous manifolds M = S 3 /T, 
we have to consider either right- or left-handed Clifford 
translations 7 6 T that act on an arbitrary unit quaternion 
q G S 3 by left-multiplication, q — > 79, respectively, right- 
multiplication, q — > q'y. The right- (left-) handed Clifford 
translations act as right- (left-) handed cork screw transla- 
tions. 

: ! ■ ■ ■ ■ ' we have carried out in lAurich et alJ |2005a) 
a systematic comparison of the CMB anisotropy with the 
WMAP data for a representative sample of homogeneous 
spherical spaces and have found agreement with the data for 
3 spaces only: T := S 3 /T* , O := S 3 /0\ and V := S 3 /I*. 



Here T* denotes the binary tetrahedral group of order 24 
generated by the right-handed Clifford translations 71 = j 
and 72 = (1 + i + j + ij)/2. O* denotes the binary octahedral 
group of order 48 generated by the right-handed Clifford 
translations 71 = (l+i)/\/2 and 72 = (1 + i + j + ij)/2. The 
last group /* corresponding to the Poincare dodecahedron 
V = S 3 /I* is the binary icosahedral group of order 120 
and is generated by the right-handed Clifford translations 
71 = j and 72 = (cr + ii + j)/2 with a = (1 + \/5)/2. In 
addition to these 3 spaces there are the 3 spaces generated by 
the analogous left-handed Clifford translations which implies 
that altogether we have to consider 6 spherical spaces. 

In order to calculate the CMB fluctuations, we have 
to expand the various contributions into the vibrations on 
the spaces M — S 3 /T which are determined by the regular 
solutions of the Helmholtz equation 



(A + = 



qeM 



satisfying the periodicity conditions 



VqeM , M lk e r 



(2) 



(3) 



Here A denotes the Laplace-Beltrami operator on S , and 
the eigenfunctions ^^''{q) are in L 2 (A4,dfi) with d/j, the 
invariant measure on S 3 . The spectrum on M is discrete, 
and the eigenvalues can be expressed in terms of the wave 



number j3 G N as E 







0* 



1 and are independent of the 
degeneracy index i = 1, . . . , r M (/3), where r M (f3) denotes 
the multiplicity of the mode (3. Only for S 3 , j3 takes all values 
in N. The allowed wave numbers j3 for the above ment ioned 
spherical manifolds are explicitly known <llkedalll995l) . see 
Table 

The relative temperature fluctuations of the CMB are 
caused by several effects which we shall compute within the 
tight-coupli ng approximation a long the lines described in 
Section 2 of lAurich et alJ i2004f) . The tight-coupling approx- 
imation leads to the Sachs- Wolfe formula Q, which incor- 
porates the various effects to be discussed below. Further- 
more, we assume here isentropic initial conditions and for 
the primordial power spectrum a scale-invariant Harrison- 
Zel'dovich spectrum 



W 2 



1 



(4) 



where a is a normalization factor which will be determined 
from the CMB data such that the angular power spectrum 
8T? fits the WMAP values in the range I = 20 to I = 45. 



3 COSMIC TOPOLOGY AND THE 

CIRCLES-IN-THE-SKY SIGNATURE 

There are infinitely many multi-connected sp herical space 
forms as described in the previous section. lAurich et alJ 
l|2005bl lar) examined the angular power spectrum ST 2 := 
1(1 + l)Ci/(2n) and the temperature two-point correlation 
function C($) for the Poincare dodecahedron T> and for 
many globally homogeneous spherical manifolds and found a 
good agreement on large scales with the WMAP data only 
for the binary tetrahedral, the binary octahedral, and the 
dodecahedral space forms. A strong suppression on large 
scales, as seen in the WMAP data, is obtained for the bi- 
nary tetrahedral space form T for fi t ot = 1.06. ..1.07, for the 
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Table 1. The eigenvalue spectrum of the considered spherical spaces S 3 and 5 3 /r. 



Ai wave number spectrum {/3} of manifold Ai multiplicity r M (/3) 
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binary octahedral space form O for fi to t = 1.03... 1.04, and 
for the dodecahedral space form T> for fitot = 1.015. ..1.02. 
This is a hint that one of these topologies could be realized 
in our Universe. Here fitot := £tot/e C rit is the energy density 
parameter, where etot denotes the total energy density of 
the universe and e C rit := |^ the critical energy density at 
the present epoch. 

In order to investigate these topologies further, we per- 
form a special search for the so-called circles-in-the-sky sig- 
nature iCornish et alJ Il993) in the CMB sky maps. This 
topological signal is caused by the multi-connectedness of 
the assumed space forms. The surface of last scattering 
(SLS) is a 2-sphere with the observer at the centre. In multi- 
connected manifolds, copies of the SLS and the observer are 
generated by the periodicity conditions If the distance 
between the observer and a copy of him ("mirror observer") 
is smaller than the diameter of the SLS, the two SLS's inter- 
sect. This intersection is a circle. Both the observer and the 
copy of the observer see this circle, but in different directions 
on the sky. Because there is really only one observer, this one 
observer should detect two circles in different directions in 
the microwave sky map. 

For a quantita tive search of these matched circles, 
ICornish et alJ (ll998T> define the quantity 



S( P ) 



(2ST a {±(P) ST b {cP + p)) 
(5Tl(J>) + 5TZ (</>)) 



(5) 



ST a (<j>) and 5Tb (4>) are the temperature fluctuations along 
two circles a and b on the microwave sky having the same 
radius. Their centres are in the directions h a and hb on the 
sky map. Both circles are parameterized by the angle <j>. 
The angle p gives the relative phase of these circles and 
( } := d(f>. By searching only for circles that are anti- 
phased, which corresponds to the minus sign in the nomina- 
tor in ©, one restricts the search to orientable topologies, 
but by searching also for circles that are phased, which cor- 
responds to the plus sign in @, one can also detect non- 
orientable topologies. The S statistic can take values in the 
interval [—1,-1-1], where S = +1 corresponds to a perfect 
agreement between the temperature fluctuations along the 
two circles and S — — 1 to completely anticorrelated temper- 
ature fluctuations. In the next section, we shall discuss the 
various effects which influence the values of the S statistic. 

In general, a complete search of matched circles requires 
a six-dimensional parameter space: four parameters for the 
directions h a and hi, of the centres of the two circles, one 
parameter for the relative phase p of the two circles, and one 
parameter for the common radius of the circles. This has to 



be done for both orientations of the circles. Such a search 
involves a huge computer time. By searching for matched 
circles of spherical manifolds, which are constructed either 
from right-handed or left-handed Clifford translations, we 
can lower the dimension of the paramet er space and conse - 
quently reduce the computational time iWeeks et al.ll2003h . 
Since we investigate in this paper three topologies which are 
globally homogeneous, we have only to search for back-to- 
back matched circles (the circle pairs are separated by 180° ) 
which reduces the parameter space to four dimensions. Fur- 
thermore, the relative phase p of the circles is given. Tables 
121 CO and 2] display the dependence of the number of paired 
circles on the radius tsls of the SLS, respectively on fitot, 
for the binary polyhedral spaces T, O and T>. In addition, 
the relative phases p and the radii of the paired circles are 
given. The column denoted by p in these Tables gives the 
relative phases of those circles which have to be added to 
the already present ones when increasing the value of fitot- 
The other circles possess the relative phases given in the pre- 
ceding row. It is possible to construct these manifolds from 
right-handed as well as from left-handed Clifford transla- 
tions leading to an identical angular power spectrum 5Tf 
and the same correlation function C(i9), but with a reversed 
sign of p in @ which alters the search for matched circles. 
Thus, one has to search for the relative phases with positive 
and negative sign. Furthermore, we can restrict our search to 
circles which are anti-phased, because the binary polyhedral 
spaces T, O and T> are orientable manifolds. 

A modificatio n of the quantity S(p), equation JSJ, 
is introduced by ICornish et alJ J2004), but which we do 
not use here. Their modified quantity, see equation (1) in 
ICornish et alJ i2004l) . is constructed by expanding the tem- 
perature fluctuations 5T a (</)) and ST/, {(f)) along the two cir- 
cles in Fourier series 8T at t,((f)) = ~%2 m T^ b exp(im^). Then 
the coefficients T^' 6 are multiplied by a factor v/|m| lead- 
ing to modified temperature fluctuations 5T a ^ (0) which in 
turn are used in ©. This procedure enhances the small 
scale structure. As will be discussed below, the available 
microwave sky maps do not have sufficient power at small 
scales, especially near to the first acoustic peak. This will 
be demonstrated by Fig.|5| where the angular power spec- 
trum 8Tf of the CMB fluctuations is compared with the 
power contained in the sky maps. With the multiplication 
by \/|m|, this deficit in power is only remedied in a statis- 
tical sense; the true cosmological signal of the SLS at small 
scales might be completely different, of course. But it is this 
very cosmological signal one needs for the detection of paired 
circles. Thus, there remains the question about how the pro- 
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cedure of lCornish et alj i 2004) affects the topological signal 
in true microwave sky maps being composed of a lot of sky 
patches. In the next section, we propose an alternative pro- 
cedure without artificially enhancing power at small scales 
but instead using a specially designed weight function. 

One can perform a very specific search by using the 
information about the group T which defines the manifold 
(S 3 /r. Fixing the radius tsls of the SLS in the case of a 
homogeneous spherical manifold, the number of paired cir- 
cles, their relative phases p and their relative orientation are 
determined. Only a common rotation of the whole configu- 
ration of the circles is left free which can be parameterised 
by three Euler angles. This fact can be used to carry out a 
simultaneous search over all matched circles by defining the 
estimator 



E(n) := sign [Y^Si • \ S -> 



(6) 



Si is the S statistic of the i-th paired circles and n is their 
total number. E varies between —1 and +1. This search 
has a four- dimensional parameter space: three parameters 
for all possible orientations of the fundamental cell and one 
parameter for the radius tsls of the SLS. The E estimator 
gives a better signal than the S statistic, if there are genuine 
matched circles; furthermore, it should be possible to detect 
with E smaller paired circles. Such a special search cannot 
be used to check all possible topologies, but if there is a hint 
towards a special topology, like the binary polyhedral spaces 
T, O and T>, it can be carried out. 



4 SEARCH FOR MATCHED CIRCLES IN THE 
WMAP SKY MAPS 

In order to perform the search for matched circles in actual 
sky maps, one has to take several points into account. The 
most important one is that one does not observe the pure 
signal of the temperature fluctuations on the SLS propor- 
tional to the scalar perturbation 4>, which displays as a scalar 
quantity a clear sign of the topology. Even ignoring the de- 
tector noise, foreground sources and all but one secondary 
contribution, i. e. not excluding the integrated Sachs- Wolfe 
effect, one has three main contributions as revealed by the 
Sachs- Wolfe formula (r(rf) := rjo — rf) 



f w - E E 

/3^3 i=l 



/ 773= 




+ 2 E E 

/3^3 i=l ''''SLS 



dr] 



''SLS 



(7) 



Here the summation runs only over those modes which ex- 
ist for a given manifold, see Tabled The first two terms in 
0, determined by the expansion coefficients $^(??) of the 
gravitational potential $ and the intrinsic perturbation in 
the radiation density p(rj), yield the effective tempera- 
ture perturbation on the SLS, i. e. the ordinary Sachs- Wolfe 



0.0 0.2 0.4 0.6 0.8 S 1.0 

Figure 1. The distribution of the S values for matched pairs 
of circles is shown for the dodecahedral space T> as the grey his- 
togram. For comparison, the distribution for S max for the simply- 
connected space <S 3 is shown as a shaded histogram, where the 
positions of the pairs are assumed to be the same as in the do- 
decahedral case. The circles have a radius of 49.8° determined by 
f^tot = 1.020. The other cosmological parameters are h = 0.70, 
O ma t = 0.28 and fi bar = 0.046. 



(SW) contribution which encodes the desired topological in- 
formation. The next term involving the spatial covariant di- 
vergence of the velocity field VlArf) is the Doppler contri- 
bution which spoils the topological signal because the pro- 
jection of the velocities in the directions of the observer and 
the "mirror" observer are not the same. The last term, the 
integral over the photon path, yields the integrated Sachs- 
Wolfe contribution (ISW). Since the photon paths towards 
the observer and the mirror observer are not identified, this 
contribution also degrades the topological signal. 

Fig. demonstrates these degrading effects. The grey 
histogram shows the distribution of S(p) (see Eq. @) with 
p — 36° for 3000 pairs of matched circles obtained from 
simulations of sky maps for the dodecahedral space T> for 
a fixed set of cosmological parameters. Without degrading 
effects, one would observe a delta- like peak at S = 1. As 
revealed by the grey histogram, the Doppler and ISW con- 
tributions lead to a broad distribution with a maximum 
around S = 0.6. The shaded histogram shows the distri- 
bution of S'max = max p {5'(p)} for 3000 pairs of circles of 
simulations of sky maps for the simply connected space S 3 . 
(For the positions of the circles we use the positions occur- 
ring in the dodecahedral simulations.) Since this space has 
no matched pairs at all, the shaded histogram provides the 
null hypothesis for the detection of a topological signature. 

Assume for a moment that the sky map has a gen- 
uine circle signature. Then one has to scan the map for 
all possible orientations of the pairs of circles. The wrong 
guesses should lead to a distribution similar to that of the 
simply-connected S 3 , whereas the few correct guesses are 
distributed as in the grey histogram. With the parameters 
chosen in Fig.0 the dodecahedral space T> has only six pairs 
(see Table Compared with the huge number of wrong 
guesses distributed as the shaded histogram, the topological 
signature can easily be washed out. 

In order to enhance the topological signal, we apply 
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Table 2. The number of paired circles for the Poincare dodecahedron T> is given in dependence on the 
radius tsls °f the SLS, respectively on the total density parameter f2tot. The parameters h = 0.70 and 
f2mat = 0.28 are held fixed. Also the relative phases p of the circle pairs are shown for the right-handed (plus 
sign) and the left-handed (minus sign) Clifford translations as well as the radii of the circles. 



T SLS 


"tot 


# of paired circles 


P 


radii of the paired circles 


< 0.33 


< 1.010 









< 0.53 


< 1.026 


6 


±36° 


< 56° 


< 0.64 


< 1.037 


16 


±60° 


< 64°, < 39° 


< 0.79 


< 1.056 


22 


±72° 


< 71°, < 55°, < 44° 


< 0.95 


< 1.080 


37 


±90° 


< 77°, < 65°, < 59°, < 38° 


< 1.05 


< 1.097 


1.3 


±108° 


< 79°, < 71°, < 65°, < 55°, < 38° 


< 1.26 


< 1.137 


53 


±120° 


< 84°, < 79°, < 76°, < 71°, < 64°, < 56° 


> 1.26 


> 1.137 


59 


±144° 


^ 84°, ^ 79°, > 76°, > 71°, ^ 64°, ^ 56°, ^ 9° 



Table 3. The same as in Table [5] for the binary octahedral space O. 



TSLS 




# of paired circles 


P 


radii of the paired circles 


< 0.40 


< 1.015 









< 0.53 


< 1.026 


3 


±45° 


< 45° 


< 0.79 


< 1.056 


7 


±60° 


< 66°, < 55° 


< 1.05 


< 1.097 


16 


±90° 


< 76°, < 70°, < 55° 


< 1.18 


< 1.121 


20 


±120° 


< 80°, < 76°, < 66°, < 44° 


> 1.18 


^ 1.121 


23 


±135° 


^ 80°, ^ 76°, > 66°, > 44°, ^ 6° 



8000 




100 200 300 400 500 I 600 

Figure 2. The angular power spectrum <5T ; 2 of the WMAP anal- 
ysis (black) of the first year data is shown . In addition the <5T 2 
spectra are shown obtained from the ILC map (red), the TdOH 
map without (blue) and with (green) a Wiener filter applied. 



a special weight function to the sky maps before carrying 
out the search for matched pairs of circles. Furthermore, we 
carry out the search for all circles of a given topology simul- 
taneously, see equation (|SJ, which should lead to a larger 
signal due to the correlation. 

Before discussing our weight function, we describe 
which angular scales seem to be preferable for such a search. 
Only maps which cover the complete sky can be used, which 
restricts us to the WMAP observations. In the following we 
concentrate on three sky m aps derived from the WMAP 
observation: the ILC map of iBennet t et alJ J2003I) and the 
two maps of iTeemark et all T2003T) with and without the 



application of a Wiener fil ter, respectively. The quality of 
these maps is discussed in jTe emark et al]l20oH : lEfstathioul 
12004 lEriksen et afll2004 120051) . In order to reveal the an- 
gular scales at which these maps can be used for a cir- 
cle search, we compare the angular power spectra STj 2 ob- 
tained from these maps using HEALPix (HEALPix web-site: 
http:/ /healpix.jpl. nasa.gov I with the angular power spectra 
obtained f rom a direc t analysis of the WMAP first year data 
(Ben nett et alJ 12003) . All four angular power spectra are 
shown in Fig. |2] Although it would be preferable to carry 
out the circle search near to the first acoustic peak, Fig. 
[5] reveals that the three derived maps do not have angular 
power spectra matching the first acoustic peak as observed 
in the WMAP spectrum. Only at scales larger than the first 
acoustic peak (I < 100) do the four spectra agree. Thus 
we cannot search for the topological signature on angular 
scales corresponding to the first acoustic peak. Instead, we 
introduce our advertised weight function which takes into 
account only the scales with I ^ 70. 

The mean value of the angular power spectrum 5Ti := 
(/(Z+l)Cy (2tt)) 1/2 is shown in Figs.|3|and|llfor the S 3 model 
and for the Poincare dodecahedral space T>, respectively. We 
also show the individual contributions to the total temper- 
ature fluctuations as discussed below equation Q , i.e. the 
ordinary Sachs- Wolfe contribution (SW) encoding the topo- 
logical signature, and the Doppler and the integrated Sachs- 
Wolfe (ISW) contributions, which both degrade the topolog- 
ical signal. The topological signal should be well preserved 
in those i-ranges, where the SW contribution dominates the 
other two, i. e. in the large scale I < 30 region and the region 
around the first acoustic peak with 100 < I < 300. However, 
as revealed by Fig. [5] neither the ILC map nor one of the 
TdOH maps have enough power up to I ~ 300 in order to 
give a reliable map for the temperature fluctuations around 
the first acoustic peak. 



© 2005 RAS, MNRAS OOO.ITHTuI 



6 R. Aurich et al. 



Table 4. The same as in Table |2] for the binary tetrahedral space T. 



TSLS 




# of paired circles 


P 


radii of the paired circles 


< 0.53 


< 1.026 









< 0.79 


< 1.056 


4 


±60° 


< 55° 


< 1.05 


< 1.097 


7 


±90° 


< 71°, < 55° 


> 1.05 


^ 1.097 


11 


±120° 


^71°,^ 55 ,>7° 
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Figure 3. The angular power spectrum <5T; is shown for an 
S 3 model with Qtot = 1 020 together with the ordinary Sachs- 
Wolfe (SW), the Doppler, and the integrated Sachs- Wolfe (ISW) 
contribution. 




2 5 10 2 5 10 2 I 

Figure 5. The weight function wi is shown for the <S 3 space (full 
curve) and for the Poincare dodecahedral space T> (dotted curve). 
In this figure, the cut-off is i cu t = 300 and not l cu t = 70 as it is 
set in the following computations. 
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Figure 4. The same quantities are shown as in Fig. El but for 
the Poincare dodecahedral space T>. 



In order to enhance the topological signal, we decom- 
pose the sky maps into the expansion coefficients ai m with 
respect to spherical harmonics Yi m (n), multiply the coef- 
ficients ai m by weight coefficients wi, and generate from 
them a new map. For the weight function u>i, we choose 
for I = 2, . . . , l cut 

wi := e X p{-(«5r max -5TD 2 /5r i L x } (8) 



with 

srr := 5T t sw - ^r ; Dopplcr + <5T/ SW ) 

and 

<5Tmax := max {ST™ ,1 = 2,..., l cut } 

For I > 1^, we set w t = 0. The coefficients <5Tf sw , 
<5y,Dop P ier and (5T/ SW are computed using the correspond- 
ing terms in equation Q. This weight function wi is unity 
for the I value with the largest SW contribution relative to 
the other two contributions with / ^ Z cut . The modes are the 
more suppressed the more important the Doppler and the 
ISW contributions are. In Fig. [5] we show an example for 
the Poincare dodecahedral space T> with the cosmological 
parameters of Fig. A significant contribution comes from 
large scale modes with I < 20. The strong decline of wi from 
I = 20 to I = 40 is mainly due to the growing Doppler term 
(see Figs. and Around I = 60. .70, the lOi's are close 
to zero because the Doppler contribution is larger than the 
SW contribution. Thus, we avoid a Gibb's like phenomenon 
by giving l cut a value just within this region. In the follow- 
ing analysis, we use l cut = 70. The contributions from the 
acoustic peak would also be important as revealed by Fig. |5] 
however, as discussed above, no reliable maps up to I ~ 300 
are at hand. 

The weight function wi differs significantly for the 
simply-connected space S 3 and a multi-connected space like 
T> at small values of I. This is due to the different eigenvalue 
spectra of these systems (see Table . The modes I = 2 
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1.010 1.015 1.020 D, tot 1.025 

Figure 6. The £ estimator is shown for two simulated sky maps 
for the Poincare dodecahedral space T> and for the simply con- 
nected <S 3 . The light grey curves show the corresponding values 
obtained from the Wiener filtered TdOH map for the left- and 
right-handed dodecahedral spaces T>. 




r 1 

1.03 1.04 1.05 O tot 

Figure 7. The £ estimator is shown obtained from simulated 
sky maps for the binary octahedral space O and for the simply 
connected S 3 . The light grey curves show the corresponding val- 
ues obtained from the Wiener filtered TdOH map for the left- 
and right-handed binary octahedral spaces O. 



and I = 3 are strongly suppressed compared to the simply- 
connected case being a welcomed property in view of the 
current discussion of the validity of the m aps with respect 
to the lowest modes iSchwarz et al.ll2004f) . For higher val- 
ues of I, the weight function wi of multi-connected spaces 
fluctuates about the values for the simply-connected case. 

In the following analysis, we use the specially designed 
S measure with respect to the space form for which the 
topological signal is looked for. This means that in the case 
of the Poincare dodecahedral space T>, the weight function 
wi is computed from the eigenvalues of T>, and all circles are 
used simultaneously occurring for that space for the cho- 
sen value of fitot. This is the best method applied up to 
now for the search of a given multi-connected space form. 
In Fig. |S| we show the results for two sky simulations for 
the dodecahedral space T> where all modes up to /3 = 185 
have been included which should suffice to obtain the contri- 
butions required by the weight function. One observes pro- 
nounced maxima near the values of fitot which are selected 
for the simulation. In addition, the result for a simulation of 
the simply-connected S 3 is shown displaying no large max- 
ima. For very small circles one can by chance find a better 
match than for larger ones. Thus, as a general trend, one 
observes a decline of E with increasing values of fi to t since 
then the radii of the circles increase. This behaviour is re- 
vealed by the S 3 simulation which has, of course, no circles 
at all. As grey curves the results belonging to the Wiener 
filtered TdOH map are shown in Fig. HJ see below. Fig. Q 
compares the E estimator for the octahedral space O with 
the one obtained from a simply-connected S 3 space. The re- 
sult for the Wiener filtered TdOH map is also shown. This 
demonstrates that our E measure yields a sufficiently large 
topological signal despite the "contaminating" Doppler and 
ISW contributions. 

In the considered interval of fitot , the Poincare dodeca- 
hedral space T> produces six matched pairs of circles on the 
sky (see Table However, with increasing values of fitot, 
the number of pairs of circles increases until the maximal 
possible number of 59 is reached. In general, the maximal 



0.90 




0.40 ^- 1 

1.05 1.06 1.07 tt tot 1.08 

Figure 8. The £ estimator is shown for two simulated sky maps 
for the binary tetrahedral space T and for the simply connected 
S . The light grey curves show the corresponding values obtained 
from the Wiener filtered TdOH map for the left- and right-handed 
binary tetrahedral spaces T. 

number of paired circles for a manifold <S 3 /T is — 1, where 
TV is the order of the group. At those values of fitot, where 
the number of pairs of circles makes a jump, the E estimator 
jumps accordingly to larger values, which is thus not a sign 
of the detection of a multi-connected space form, but simply 
an artefact of the increasing contributing number of pairs to 
the E estimator. This is the case for the binary tetrahedral 
space T at fitot = 1.056 where the number of circle pairs 
jumps from 4 to 7 (see Table [IJ. In Fig. [8] the E estimator 
is shown for two sky simulations for the binary tetrahedral 
space T having a value of fitot = 1.054 and fitot = 1.065, re- 
spectively. Both show a maximum near to the expected value 
of fitot, but in addition a superimposed peak at fitot = 1.056 
is visible where the number of circle pairs jumps from 4 to 
7. The same peak is seen in the S 3 simulation which has 



© 2005 RAS, MNRAS OOO.HHTuI 



8 R. Aurich et al. 



no matched pairs of circles. As in Fig. |S| the grey curves 
display the results belonging to the Wiener filtered TdOH 
map. (The peak at fi to t = 1.056 is also seen in Figs. I13| and 

M) 

Now, we apply the E estimator to the Poincare dodec- 
ahedral space T>, the binary octahedral space O, and the 
binary tetrahedral space T using thre e sky maps obt ained 
from the WMAP measurements feennett et alJ l2003). The 
maps are constructed using different combinations of the 
data obtained in different frequency bands in order to sup- 
press foreground sources. Different combinations are used for 
different sky patches. This is not ideal for the search of pairs 
of circles, since this leads to an inhomogeneous structure in 
the maps. The E estimator correlates the values on the cir- 
cles which are then obtained by different combinations of the 
original data. The WMAP team explicitly states that their 
ILC map should not be used for cosmological studies. But 
due to a lack of an alternative, we apply, with this drawback 
in mind, the S estimator to the IL C map produced by the 
WMAP team ([Bennett et al] |2003'l. and to the two maps of 
iTegmark et al l <2003ft ~ 

In Figs. [5] and 1101 the results for the E estimator are 
shown for the Poincare dodecahedral space D for the right- 
and left-handed space form, respectively. For the dodecahe- 
dral space T>, the best match to the angular power spectrum 
measured by WMAP is obtained around fitot = 1.018 in 
lAurich et all J2005al) . Is there something special in Figs. |U] 
and 1101 near to this value of fitot? Whereas the left-handed 
variant shown in Fig. If 01 displays only the smooth decline 
of E due to the increasing size of the circles with increas- 
ing values of fitot, the right-handed space form shown in 
Fig. [5] displays a clear maximum around fitot = 1.015 in 
all 3 maps which is par ticularly distinct in the two maps 
of ITegmark et all (|2003). One is tempted to interpret the 
peak at fitot = 1.015 as an indication that the shape of the 
universe is that of the right-handed Poincare dodecahedron 
T). However, the maximum seen in Fig. [0] is far less pro- 
nounced than the maxima obtained from our simulations 
of sky maps for the dodecahedral space T>, as can be seen 
by a comparison with Fig. |H|where the simulations (full and 
dotted curves) are presented for two different fi to t-values to- 
gether with the Wiener filtered TdOH values (grey curves). 
There, one observes that the height of the maximum is of the 
same order as the amplitude of the S 3 simulation possessing 
no matched circle pairs. However, these simulations are not 
based on a patchwork of sky maps, and it is expected that 
such a patchwork degrades the signal for a multi-connected 
space. But since its is unknown by how much the signal is 
degraded, we cannot draw firm conclusions from this struc- 
ture at this moment due to the limited quality of present 
CMB data. 

A possib le hint for the left-han ded dodecahedral space 
T> is found bv lRoukema et al] (120041) who find a slight signal 
for six circle pairs with 11 ± 1° in the ILC map which is not 
revealed distinctly in Fig. 1101 

The E estimator is shown for the right- and left-handed 
binary octahedral spaces O in Figs. Ill) and ll2l respectively. 
For this space form, the best match to the WMAP angular 
power spectrum is o btained for values around fitot = 1.038 
\ Aurich et afll2005alV The right-handed space displays only 
the decrease due to the increasing size of the circles. The 
left-handed case displays a minimum near fitot = 1-04 and 
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Figure 9. The results for the £ estimator are shown evaluated 
for the right-handed Poin care dodecahedral spa ce T> for the ILC 
map and the two maps of Tcgmark et al. (2003). 
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Figure 10. The same as in Fig. [5] is shown for the left-handed 
Poincare dodecahedral space T>. 



E increases for higher values of fitot only to values typical 
of the S 3 simulation. 

Figs. ll3l and ll4l show the E estimator for the tetrahedral 
space T again for the right- and th e left-handed variant, 
respectively. In lAurich et alJ (l2005al) it is found that the 
suppression of large scale power agrees with the WMAP 
measurement for fitot in the range 1.06 . . . 1.07. The large 
peak observed around fitot = 1.056 is caused by the increase 
of the number of circle pairs which jumps from 4 to 7 (see 
Table EJ. The right-handed space form displays large values 
around fitot — 1.068. As in the case of the dodecahedron, 
the maximum is far less pronounced than in the simulation 
of a T sky as revealed by Fig. El an d the remarks concerning 
the dodecahedron's maximum applies also to this case. 



5 CONCLUSION 

We performed a thorough search for a circle-in-the-sky sig- 
nature in CMB maps under the hypothesis that the topol- 
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Figure 11. The same as in Fig. [5] is shown for the right-handed 
binary octahedral space O. 
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Figure 14. The same as in Fig. [5] is shown for the left-handed 
binary tetrahedral space T . 
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Figure 12. The same as in Fig. |^]is shown for the left-handed 
binary octahedral space O. 
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Figure 13. The same as in Fig. O is shown for the right-handed 
binary tetrahedral space T. 



ogy of the universe is given by one of the 3 spherical spaces 
known as 23, T, and O. To this purpose, we introduced the 
£ estimator © which is applied to filtered sky maps us- 
ing the weight function JHJ in order to enhance the ordinary 
Sachs-Wolfe contribution. A possible signal was found in 
the case of the right-handed Poincare dodecahedral space 23 
for fitot — 1.015, see Fig. |U] which can be interpreted as a 
hint that our universe possesses the topology of 23. A simi- 
lar possible signal occurred in the case of the right-handed 
binary tetrahedral space around fitot — 1-068, see Fig. 1131 
However, one must keep in mind that all three sky maps, 
derived from the WMAP measurements in five frequency 
bands, are a patchwork and thus possess a complicated sig- 
nal to noise behaviour. Since the use of these incomplete 
sky maps certainly degrades a possible topological signal, 
and in order to avoid a "false negative", we can only re- 
port our findings without drawing a final conclusion about 
the true topology of the Universe. Furthermore, we believe 
that lCornish et alJ 112004!) overstate their result of not find- 
ing paired circles since they also rely on the same CMB data. 
The question whether our Universe possesses a non-trivial 
topology is to our opinion not yet decided. 
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